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Abstract — In this paper, by assuming equiprobable centro- 
symmetric constellation adopted for modulation, the mathe- 
matic expectation of the sample auto-correlation matrix obtained 
through the least-squared channel estimation is derived under 
a rigorous model of the doubly selective fading channel for 
orthogonal frequency division multiplexing (OFDM) systems. 
Based on the derived expectation, the bias-property of the sample 
auto-correlation matrix is analyzed with respect to the total 
squared error (TSE) and the eigen-structure. And according to 
the analyses, when the sample auto-correlation matrix is used 
to estimate the signal and noise subspaces, BPSK may not be 
a good choice for modulation of pilot tones, which contradicts 
with the experience of most current OFDM systems. Numerical 
results demonstrate this bias-property and the analyses on it. 

Index Terms — Bias, Sample auto-correlation matrix. Doubly 
selective fading channels, OFDM, Total squared error, Eigen- 
structure. 

I. Introduction 

The auto-correlation matrix of the channel frequency re- 
sponse (CFR) plays a very important role in the channel 
estimation for orthogonal frequency division multiplexing 
(OFDM) systems. Many channel estimation algorithms are 
based on this statistics, for example, the linear minimum mean- 
squared error (LMMSE) estimator and its optimal low-rank 
approximations [1], the MMSE estimator exploring both time 
and frequency correlations [2], the two-dimensional Wiener 
filtering [3], and those algorithms based on parametric channel 
modeUng [4] [5] [6]. 

However, the true auto-correlation matrix of CFR is un- 
known in real applications, hence, the sample auto-correlation 
matrix is used in stead and usually obtained through the least 
squared (LS) channel estimation. For fixed or slowly moving 
radio channels whose Doppler spreads are relatively small, the 
channels reveal a feature of block-fading [7], which means 
the channels are static or nearly static intra-symbol. Under the 
circumstances, the sample auto-correlation matrix is unbiased 
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if noise has not been involved. However, for fast moving radio 
channels, intra-symbol fading becomes so prominent that inter- 
carrier interference (ICI) takes effect by not only degrading 
the performance of OFDM systems [8], but also biasing the 
sample auto-correlation matrix. In [9], the bounds of the ICI 
power is derived. 

In this paper, we focus on the bias-property of the sample 
auto-correlation matrix for OFDM systems in doubly selective 
fading channels. Based on an accurate channel modeling, 
we derive the mathematic expectation of the sample auto- 
correlation matrix by assuming equiprobable centro-symmetric 
constellations adopted for modulation. Further analyses, in- 
cluding the total squared error (TSE) and eigen-structure, are 
conducted thereafter 

This paper is organized as follows. In Section|IIl the OFDM 
system and channel model are introduced. Then, derivations 
of the bias of the sample auto-correlation matrix and related 
analyses are presented in Section |lll] Numerical results appear 
in Section |IV] Finally, Section [V] concludes the paper. 

Notation: Lowercase and uppercase boldface letters denote 
column vectors and matrices, respectively. 0, (•)*, (•)^, and 
II • IIj;- denote Hadamard product, conjugate, conjugate trans- 
position, and Frobenius norm, respectively. E{-) and Ex{-) 
represent the expectation and conditional expectation with 
respect to x, respectively. [A]i,j and [a]i denotes the (i,j)-th 
element of A and the i-th element of a, respectively. diag{a) 
is a diagonal matrix by placing a on the diagonal. 

II. System Model 

Consider an OFDM system with a bandwidth of BW = 
1/T Hz (T is the sampling period). N denotes the total number 
of tones, and a cyclic prefix (CP) of length Lcp is inserted 
before each symbol to eliminate inter-block interference. Thus 
the whole symbol duration is Tg = {N + Lcp)T. 

The complex baseband model of a linear time-variant mo- 
bile channel with L paths can be described by [10] 



hit,T)=J2hiit)Sir-TiT) 



(1) 



1=1 



where r; e 7?^ is the normalized non-sample-spaced delay 
of the ^-th path, and hi{t) is the corresponding complex 
amplitude. According to the wide-sense stationary uncorre- 
cted scattering (WSSUS) assumption, /i;(t)'s are modeled as 
uncorrected narrowband complex Gaussian processes. 

Furthermore, by assuming the uniform scattering environ- 
ment introduced by Clarke [11], /i;(t)'s have the identical 
normalized time correlation function (TCF) for all /'s, thus 
the TCF of the Ts path is 

n,i(At) =afJo(2^/dAt) (2) 

where af is the power of the l-th path, fd is the maximum 
Doppler spread, and Jo(-) is the zeroth order Bessel function 
of the first kind. Additionally we assume the power of channel 
is normalized, i.e., '^'i ~ ^■ 

Assuming a sufficient CP, i.e., Lcp > L, the discrete signal 
model in the frequency domain is written as 

y/ (") = H/ (n)x/ (n) + Uf (n) (3) 

where x/(n), y/(n), n/(n) G C^^^ are the n-th transmit- 
ted and received signal and additive white Gaussian noise 
(AWGN) vectors, respectively, and H/(m) e C^^^ is the 
channel transfer matrix with the {k + k)-tii element as 

N-l L 
m=0 1=1 

where hi{n,m) — hi{nTs + {Lcp + m)T) is the sampled 
complex amplitude of the l-t\\ path, k and v denote frequency 
and Doppler spread, respectively. Apparently, as H/(n) is 
non-diagonal, ICI is present. In fact, when fdTs < 0.1, the 
signal-to-interference ratio (SIR) is over 17.8 dB [12]. 

III. Analysis of the Sample Auto-Correlation 
Matrix of the Channel Frequency Response 

Usually the auto-correlation matrix of the channel frequency 
response is obtained through the LS channel estimation. For 
systems adopting pilot-symbol assisted modulation (PSAM) 
[1], only pilot symbols, denoted as yp{n) e C^^^, are used to 
perform LS channel estimation. For systems employing data- 
aided channel estimation [13], all tones can be utilized to assist 
channel estimation if they have been demodulated. Then, 

hpjs{n) = Xp^(n)yp(n) 

= X~\n)Upin)Xp{n) + X-\n)np{n) (5) 

where Xp(n) = diag{xp{n)) is a diagonal matrix consisting 
of pilot or demodulated symbols, and the noise term is 
np{n) ^C^{0,ctIIn). 

The auto-correlation matrix of CFR referred in most litera- 
tures (see citations in section IJl is assumed of the form as 

Rp = Fp,,DpF^, (6) 

where Fp i- e C^^^ is the unbalanced Fourier transform 
matrix, defined as [Fp^r]kj = e"-'^'^'^^'/^. Dp e C^^^ is a 
diagonal matrix with the l-th diagonal element [Dp]ij = erf. In 



TABLE I 

Values of p for Typical Modulations 



Type of Modulation 


BPSK 


QPSK 


16QAM 


64QAM 


P 


1 


1 


L889 


2.685 



fact. Dp is the auto-correlation matrix of the channel impulse 
response (CIR). 

With the LS estimated CFR h.p,is{n), the sample auto- 
correlation matrix of CFR is constructed in most cases as 

1 ^* 

* n=l 

In most literatures, (|7]l is treated as an unbiased estimator 
of ^ when noise in not involved in the channel. However, 
substituting (|5]l into d?), we have dHJ, shown at the bottom of 
the next page, where Nt is the total number of observations. 
Therefore, its mathematic expectation is 

1 ^' 

EiRpj,) = E E[KiAn)h^js(.n)] (9) 

* Tl=l 

By assuming that channel fading is independent of thermal 
noise, the additive term in (|9]l is 

E[hp.UnKAsH] = 

E[X- 1 (n)Hp (n)xp (n)x^ (n)Hf (n)X;^ (n)] 
+alE[X-\n)X-"{n)] (10) 

In addition, we assume that the elements of Xp(n), n — 
1, . . . , Nt, are independent of each others, which is usually 
the case for pseudo-random pilot sequences or modulated data 
symbols. And so are the elements of Xp{n). Besides, we also 
assume a normalized centro-symmetric constellation set with 
equiprobability. Therefore, 

E{[^p{n)]k) = (11) 

£;([Xp(ni)]fcJx;(n2)]feJ = - n2)5(fci - fca) (12) 

i?([xp(ni)]-^i[x;(n2)]fc,') = pS{ni - n2)S{k, - k2)(l3) 

where p is a constant related to the type of modulation. For 
MPSK, p — I, and for MQAM, p depends on the size of 
constellation, i.e., 

^ 3M ^ ^ 

fcl=0 fe2=0 

[{2ki + 1 - VM)2 + (2fc2 + 1 - VM)^]-! (14) 

where M = 2^™, m =1,2,.... Values of p for several typical 
modulations are evaluated in Table H] 

For centro-symmetric constellations, according to ( fT2b and 
(O, the first term on the right-hand side of ( fTOl i is expressed 
as follows through some manipulations (see Appendix lAb. 

i?[X;i(n)Hp(n)xp(n)xf(n)Hf(n)Xp-^(n)] = 

T Fp,,DpF^, + (1 - ^o)ptIIp, BPSK 

^oFp,,DpF^, + (1 - ^o)palIp, others 



where — J^iLo '^f — 1' ^^'^ is a scalar, defined as 

N-lN-l 

^0 = ]^ E E M^nfdim - q)T) (16) 

m=0 q=0 

and T is a Toeplitz matrix, defined as 

[^]ki,k2 = E E M'^T^fd{mi~m2)T) 

mi— 7712— 

xe^-'^'^'^'^^^'^^)'^™^^'"^-' (17) 

Substituting ( fT3] l into dTol l. the second term on the right- 
hand side of ( [Tol l is expressed as 

alE[lL-\n)y.-" in)] = pa^Ijv (18) 

With ([l5]l([T8]l([T0ll, (ID is rewritten into 

^ r T © Fp,,DpF^, + [(1 - eo)7 + llpfT^Jjv, ^^^^.jg. 

where 7 = -j- is the signal-to-noise ratio (SNR). 

If considermg the noise term presented in the LS estima- 
tion, the expected auto-correlation matrix, i.e., should be 
modified to 

HpJs = Fp,,DpF^^ + pgIIn (20) 

Compared with ( |20l ), ( fT9] l is no longer an unbiased estimator, 
and for most types of modulations except BPSK, its bias is 

AR — Rpjs — E{Rp^is) 

= (1 - eo)Fp,.DpF^, - (1 - ^o)p4'^n (21) 

Hence, the corresponding TSE is defined as 

L 

TSE(Rp,i,) ^ ||AR|||, = (1 - Co)' Y.^h - polf (22) 

where Aj's are the eigenvalues of Rp. 

Although dTji reveals the property of bias, according to ( |T9] l. 
using £'(Rp is) to estimate the signal and noise subspaces is 
immune when BPSK is not adopted. In fact, let the eigenvalue 
decomposition (EVD) of Rp,;s is [14] 

Rp,,, = Fp,,DpF^, + pallN = UrU^ 

where U is unitary, and F is diagonal with \r]k.k — A; +pcr,^j. 
Then, the EVD of £'(Rp ^s) can be expressed as 

E(%,is) = eoFp,,DpF^, + [(l-eo)7+l]P'^'liv = Uf 

where [f]fc,fe = ?oA, + [(l-eo)7+l]/"^'. and [r]fe,fe- [f ]fc,fc ^ 
(1 — ^o)(A, — pcrfj)- Therefore, only the eigenvalues are under- 
mined by the interferences introduced by the doubly selective 
channels, leaving the corresponding eigenvectors unimpaired. 




y{dB) 



Fig. 1. Evaluating r versus f^Ts and 7. 

Nevertheless, if BPSK is adopted, not only the eigenvalues 
are undermined by the interference, but also the eigenvectors, 
of which the signal and noise subspaces consist, are deterio- 
rated due to the presence of T. 

Besides, from ( fT9t , we notice that the ICI caused by the 
doubly fading channels boosts the relative level of noise at a 
ratio of r = ^q"^ + (^q"^ — 1)7 = (1 + 7)'Co''^ ~ 7- Since when 
jdXs < 0. 10, Co is a strictly monotonically decreasing function 
with respect to /d, r is strictly monotonically increasing with 
respect to /d, meanwhile, as Co<l^ ^ is strictly monotonically 
increasing with respect to 7. Fig[T| plots r versus jdTs and 7, 
where /^T^ is the normalized Doppler. 

IV. Numerical Results 

The OFDM system in simulations is of BW = 10 MHz 
(T = 1/BW = 100 ns), N = 1024, and L^p = 128. 
Two 3GPP E-UTRA channel models are adopted: Extended 
Vehicular A model (EVA) and Extended Typical Urban model 
(ETU) [15]. The excess tap delay of EVA is [0, 30, 150, 310, 
370, 710, 1090, 1730, 2510] ns, and its relative power is 
[0.0, -1.5, -1.4, -3.6, -0.6, -9.1, -7.0, -12.0, -16.9] 
dB. For ETU, they ai-e [0, 50, 120, 200, 230, 500, 1600, 
2300, 5000] ns and [-1.0, -1.0, -1.0, 0.0, 0.0, 0.0, -3.0, 
—5.0, —7.0] dB, respectively. The classic Doppler spectrum, 
i.e., Jakes' spectrum [10], is applied to generate the Rayleigh 
fading channel. 

In Fig 121 Nt = 10^ OFDM symbols are collected to 
construct the sample auto-correlation matrix to approach its 
expectation. Three different types of modulation, including 
QPSK, 16QAM and 64QAM, are considered as well as a wide 
range of SNR's and Doppler's. It is shown that the numerical 
results coincide with the analytic results ( |22] | rather well. 
According to the figure, TSE's for higher Doppler's are much 

*Tliis condition is implicit for applicable OFDM systems to maintain the 
power of ICI within a tolerable range [12]. 



1 ^* 



(8) 
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Fig. 2. Evaluating TSE of the sample auto-correlation matrix for different 
modulations, Doppler's and SNR's when Nt = 10®. 



larger than lower ones, which follows ( |22] | that ||AR|||, is 
strictly monotonically increasing with respect to fd- Moreover 
since 



ailARii. =-2(l-eo)'(^"Ef=i'^f-ip) <0 for large 



dp 

N, TSE is decreasing with respect to p, which is demonstrated 
by the figure, too. Besides, for a given channel, the effect 
of different modulations on TSE is not as significant as of 
Doppler's. 



V. Conclusion 

In this paper, the bias-property of the sample auto- 
correlation matrix for OFDM systems in doubly selective 
fading channels is derived and analyzed. Through the analyses, 
it is found that the TSE of the sample auto-correlation matrix 
is increasing with respect to the Doppler but decreasing with 
respect to the order of modulation. Although the intra-symbol 
fading biases the expectation of the sample auto-correlation 
matrix, the eigenvectors are not affected at all when BPSK is 
not adopted, which allows subspace-based channel estimation 
algorithms. However, when BPSK is adopted, like most current 
OFDM systems did, not only eigenvalues but also eigenvectors 
are undermined by ICI due to the intra-symbol fading. 

Appendix A 
Derivation of (fTsI ) 

Let A = Xp i(n)Hp(n)xp(n)x^(n)H^(n)Xp^(n), then 

N N 
qi = l q2 = l 



mi=0 h = l 
N-1 L 



;^ E E K{n,m2)t 



-]2TT{{k2-q2)m2+k2Ti^)/N 



m2=0 l2 = l 



iV — i > 



Therefore 

{N N 
E E [a;-^(n)a;^;(n)a;5,(n)a;;jn)] 
gi = l 92=1 

(N-1 L N 
^ E E^'i("'"^i)^"''""''"''^"'^''''^^^'^ 
mi=0/i=l / 

N-1 L 

\ m,2=0l2 = l / ) 

With (fTTli-(fT3]). we evaluate Ex[-] under the following cases. 

1) When ki = k2 ~ qi — q2 — k, 

Ex [xl^ {n)x~^* {n)xk{n)xl{n)\ = 1 

2) When ki = k2 — k, qi — q2 — q, ky^q. 

Ex [x'^^{n)x'^*{n)xq{n)xl{n)] = p 

3) When fci = qi, A;2 = 92, ^I7^fc2, 

[a;fclH?^)a;fe,*Na;fel(")42(")] = 1 

4) When fci = 92, = 91,^17^^2, 

' 1, BPSK 
0, others 



Ex [x^l{n)x^*{n)xk2{n)xl^{n)\ 
5) Others 



Ex [x^l{n)x,^*{n)xq^(n)x*q^{n)\ = 



Note: To proved, let Xki{n) = re^^, then x^^ {n)x*f._^{n) = 
g-i26i equiprobable centro-symmetric constellations, 9 

and —6 are equiprobable, so i?e(e~^^*) = for all centro- 
symmetric constellations, e.g., MPSK (M ^ 2) and MQAM. 
Only for BPSK, Ee{e-^'^% = 1. 

According to [U and |2]i, ( l23T l. shown at the bottom of this 
page, is obtained. Further, notice that 

Eh[hi^{n,'mi)h*i^(n,m2)\ = (Tf^5{li-l2)Jo{2Tifd{mi-m2)T) 
(l23T l is expressed as 

/ L \ / N \ 



-m2)/N 



^1 = 1 



(N-1 N-1 \ 
]^ E E ^o(2vr/,(mi-m2)T)e-^"2-'=(™^-™^)/^ 
mi— 0m2— / 

^(l-/')(E^n (4 E E ^o(2^/d(mi-m2)r)] 



mi— 7712—0 



Since J2q=i e^'2^«(™i-™2)/^ = iV(5(mi - ma), we have 

E {[A]k^k) = P'jI + (l - pWUo 
where cr^ — X^^Li — 1' is defined as 

^0 " 1^ E E ^o(27r/rf(mi - m2)T) 



(25) 



mi— m2— 



Similarly, according to [3]l-|5]l, (l24l) is obtained and shown at 
bottom of this page. Therefore, 



1=1 

L 



others 
(26) 



where fci^fca, and 

JV-l w-i 

-EE ^o(27r/d(mi-TO2)T)e-^"2-('=i-'=^)('"i+'"^) 



Consequently, with ( IZSl l and ( |26] l, we have 

[X; 1 (n)Hp (n)xp (n)x^ (n)H^ (n)X;^ (n)] = 
j Y Fp,,DpF^, + (1 - eo)p'T2lp, 
I ^oFp.rDpF^, + (1 - eo)p^T,'lp, of/ieri 

where Y is a Toeplitz matrix, defined as ['^]ki,k2 — '^ki-k2- 
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mi— m2— 



N N-1 L N-1 L 

q=l mi=0 h = l ni2=0 l2=l 
N-1 L N-1 L 

E E E E^'^['^'^(-'™i)^'*.("'-2)]^"''''^"^"^""=^'^'' (23) 



mi— /i — 1 m2— /2 — 1 



N-1 L N-1 L 

1^ E E E E £^'i[/»ii(">"ii);i,;(n,m2)]e-J2'^[('=i-'=2)('"i+"2)+'^i^'i-'=2^'2l/^, 

TP f\ \] ^ / mi=0 Zi = l m2=0 ;2 = 1 

^ \.l-^\ki,k2) — \ AT-i L N-1 L 

172 E E E E £;,,[/i,i(n,mi)/i,;(n,m2)]e-J"2-('^i-'i-^-^-'2)/^, others 



mi— li — l m2— /2— 1 



(24) 



